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ABSTRACT 

There is both theoretical and empirical evidence that the initial mass function (IMF) 
may be a function of the local star formation conditions. In particular, the IMF is 
predicted to flatten with increasing local luminosity density pi, with the formation of 
massive stars being preferentially enhanced in brighter regions. In R136, the bright 
stellar cluster in 30 Doradus, the IMF gradient is dV/dlogpi = 0.28 ± 0.06, where F 
is the slope of the IMF. If such IMF gradients are indeed general features of galaxies, 
this implies that several previous astrophysical measurements, such as the surface mass 
densities of spirals (obtained assuming constant mass to light ratios), were plagued 
by substantial systematic errors. In this Letter, calculations which account for possible 
IMF gradients are presented of surface densities of spiral galaxies. Compared to previous 
estimates, the mass surface densities corrected for IMF gradients are higher in the outer 
regions of the disks. For a model based on the Milky Way but with an IMF scaled 
according to R136, the rotation curve without the traditional dark halo component 
falls with Galactocentric radius, though slower than it would without IMF gradients. 
For a second model of the Milky Way in which the IMF gradient is increased to 0.42, 
the rotation curve is approximately flat in the outer disk, with a rotational velocity 
below ~ 220 km s~^ only before the traditional dark halo component is added. For 
a third model in which substantial arm/inter arm density contrasts are additionally 
assumed, the solar vicinity mass density drops to O.IOMqpc"^, which is consistent with 
observations. These results, if generalizable to other galaxies, not only call into question 
the assertion that dark matter halos are compatible with the flat rotation curves of spiral 
galaxies, but also may clarify our understanding of a wide variety of other astrophysical 
phenomena such as the G-dwarf problem, metallicity gradients, and the Tully-Fisher 
relation. 



Subject headings: galaxies: luminosity function, mass function — galaxies: kinematics 
and dynamics — dark matter — galaxies: halos — Magellanic Clouds — galaxies: 
evolution 



^Postal address: NASA/GSFC, Laboratory for High Energy Astrophysics, Code 661, Greenbelt MD 20771, USA 
^Department of Physics, University of Maryland, College Park, MD 20742, USA 



- 2 - 



1. INTRODUCTION 

In a recent paper, Padoan, Nordlund, k. Jones (1997) claimed on theoretical grounds that 
the initial mass function (IMF) should be a function of the local temperature T of the original 
molecular clouds. Padoan et al. (1997) argued that dense star forming regions, such as those in 
starburst galaxies, should be warmer than sparser star forming regions. In fact, if the temperature 
dependence of the clouds is not drastically different from that of a blackbody, then T oc p\^^ , where 
p\ is the local mean luminosity once star formation has already started. Padoan et al. claimed that 
starburst regions should therefore have a flatter IMF and be more "top heavy." Similar reasoning 
would imply that the IMF in cooler regions of galaxies should favor low mass star formation and 
be steeper. 

In support of their star formation model, Padoan et al. (1997) noted that for T ;^ 60 K, their 
models predict a top heavy IMF similar to that found in the center of R136 (Malumuth & Heap 
1994, Brandl et al. 1996), the bright stellar cluster in 30 Doradus. Due to its proximity and the 

fact that it is the most massive H II region in the Local Group, 30 Doradus is perhaps the best star 
formation "laboratory" accessible to us. However, the relaxation time in R136 may be less than its 
age (Campbell et al. 1992), so dynamic friction may also contribute toward the R136 present-day 
mass function gradient. 

Fortunately, many other avenues of testing Padoan et al.'s model exist. The 0-star catalog of 
Garmany, Conti, &; Ghiosi (1982) shows a flattening of the IMF slope toward the Galactic center 
(cf Humphreys & McElroy 1984). This data supports Padoan et al.'s model since higher surface 
brightness regions would, on the average, yield higher temperatures and flatter IMFs. Models which 
attempt to explain correlations between local surface brightness, color, line ratios, metallicity, and 
the star formation rate have assumed luminosity-dependent IMFs (e.g., Edmunds & Phillipps 1989; 
Phillipps, Edmunds, &: Davies 1990). Several evolutionary models of inner regions of starburst 
galaxies assume low mass cutoffs or top heavy IMFs (e.g., Rieke et al. 1980; Augarde & Lequeux 
1985; Doane k, Mathews 1993; Doyon, Joseph, k, Wright 1994). Finally, independent theoretical 
arguments supporting IMF gradients range from models which are consistent with the simple form 
of the Jeans expression for the typical stellar mass in solar units of < m >oc T^/^ (e.g.; Larson 
1982; Bodcnhcimer, Tohline, k Black 1980) to much more complicated models, such as the outflow- 
regulated model of Adams k Fatuzzo (1996), which predicts < m >(x T", where 1 < a < 3/2. 

If IMFs are actually a function of pi or T, there would several important astrophysical conse- 
quences. For instance, there would be a position-dependence in the mean mass to light ratio. This 
is due to the strong dependence of the mass to light ratio upon the IMF. In R136, this makes the 
mass density function pm much different from pi (Malumuth k Heap 1994, Brandl et al. 1996) and 
complicates estimates of the total mass. Padoan et al.'s results indicate that similar effects might 
occur in spiral galaxies. If the luminosity of a star is taken as L ~ LqUi'^, where y ~ 3.5, the Jeans 
expression above would suggest the crude relation < m >oc pj*^^ and yield pm oc pj^^''^ ^''^^ — p^'^^. 
Unfortunately, previous works have assumed that IMFs are independent of time and position with. 
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specifically, oc pj-^ throughout a given spiral galaxy (e.g., van Albada et al. 1985). In this Letter, 
surface mass densities of spiral galaxies are computed, for the first time, by explicitly accounting 
for the possible types of IMF gradients that might exist if theories like those of Padoan et al. are 
correct. 

2. AN EMPIRICAL ESTIMATE OF THE R136 IMF GRADIENT 

Since position-dependent measurements in R136 of both pi and the IMF slope T (where 
dN/dM oc mF'^ is the number of stars per unit mass in solar units) have already been made, 
computing the dependence of the R136 IMF upon the local luminosity is straightforward. Doing 
this will provide a Tiscful starting point in obtaining a crude yet quantitative estimate of the possible 
types of IMF gradients that might generally exist in all galaxies including the Milky Way. 

Table 1 summarizes Brandl et al.'s (1996) results for the IMF based upon high resolution 
5-color photometry of the stars in R136 estimated to be between 2.5 and 3.5 Myrs years old. The 
right-most entry of Table 1 shows the results of performing the coordinate transformation between 
R and pi using Figure 15 of Hunter et al. (1996). Though Brandl et al. (1996) did not make 
explicit measurements of the upper and lower stellar mass cutoffs mi and to the power-law 
approximation of the IMF, Table 1 includes estimates of their dependences upon the local surface 
brightness. The lower mass limits were obtained from the peaks of Brandl et al.'s mass functions, 
while the upper limits were taken from the highest masses observed per radius bin. Both log^g(mi) 
and \ogiQ{m^) are found to decrease by ~ 0.2 with each successive increase in radius. Brandl et al. 
(1996) performed completeness corrections, so the depletion of low mass stars in all but the outer 
regions of R136 is presumably real. The results of performing a linear fit of the IMF parameters 
of R136 to logio(/9i) are shown in columns 2-7 of Table 2 as Model A. Uncertainties of parameters 
calculated from more than two radius bins are shown in parenthesis. 

3. DYNAMICAL PROPERTIES OF SPIRAL GALAXIES WITH IMF 

GRADIENTS 

The IMF gradient of Model A implies a surface mass density that is different from what would 
be obtained were the mass to light ratio constant. The surface mass density for Model A, if scaled 
according to the surface luminosity function suspected for the Galaxy, is shown in the top panel 
of Figure 1. The disk scale length Rq = 4.5 kpc and solar Galactocentric radius Rq = 7.8 kpc 
were taken from Kuijken &; Gilmore's (1989a) model of the Galaxy. For simplicity, pi at a given 
radius was assumed to be constant throughout a disk thickness of 575 pc. The surface brightness 
was normalized to be 22.5 LqPC^^ at R = Rq, which results in piQ = pi{Rq) = 0.037 LqPC~^. 
For each radius bin, the IMF was obtained from pi and the coefficients shown in Table 2 for Model 
A. This IMF was converted to present day mass and luminosity functions by assuming (purely 
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for simplicity) a constant star formation rate for the past 1.0 x 10^° yrs. The main sequence 

Ufetimc-himinosity-mass relationships used to obtain the mass to light ratio as a function of the 
IMF were obtained from logarithmic-linear interpolation of m > 0.8 models published by Schaller 
et al. (1992) for z = 0.02, overshooting of the m > 1.5 stars, and standard mass loss rates. For 
m < 0.8, LU=o.25 = 7.8 x lO'^L©, L |m=o.08 = 6.55 X 10 ^-^0, and L\m<o.o7 = 5.0 x 10 ^'^Lq were 
assumed. 

The surface densities both of Model A and of the constant mass to light ratio model fall off 
exponentially with increasing radius. The effective scale length of Model A is ~ 7.5 kpc, which is 
2± 1.7 times larger than that of the surface brightness function. This increase in the scale length is 
a result of the fraction of low mass stars (and the mass to light ratio) increasing with radius. 

From the surface density, other dynamical properties of the galaxy can also be calculated. The 
circular velocity (i.e., the rotation curve) corresponding to the surface density of Model A is shown 
in the second panel of Figure 1. The parameters for the bulge, spheroid, and halo were taken from 
Table 1 and Figure 5 of Kuijken & Gilmore's (1989a) model of the Galax;y. To avoid a divergent 
and unphysical total mass, the additional assumption that all components of the Galaxy terminate 
at an arbitrarily-selected maximum radius of 35.0 kpc was also made. For Model A, this results 
in a total mass of the halo, bulge/spheroid, and disk of, respectively, 2.8 x 10^"'^Mq,3.5 x W^^Mq, 
and 6.3 x W^^Mq. In comparison, the integrated disk mass of the 7v = 2.0Mq/Lq model is only 
3.2 x IO^^Mq and increases much faster with radius. For simplicity, the surface mass density of 
stellar remnants and gas was assumed throughout the disk to be 1/3 that of the stars. Because 
the halo dominates the mass distribution, the circular velocity curve (solid line) is nearly flat. 
Without the halo, the circular velocity curve falls from 185 km s~^ at i? = 2.0 kpc to 124 km s~^ 
at ii = 34 kpc. Though the surface density of Model A corrected for IMF gradients is different 
from that previously obtained for spiral galaxies. Figure 1 shows that the change is not enough to 
dramatically affect the dynamical properties of the disk, such as the circular velocity curve. 

The IMF of Model A is very negative at all radii, with r|R=i.o kpc = —'^■'2,^\r=Rq = —3.4, 
and r|j:j=i5.o kpc = —3.6. This occurs even though R136's spatially-averaged IMF is typical and its 
IMF gradient is small only because it has a luminosity density that is 10^ — 10^ times higher 
than typical regions of spiral galaxies. For comparison. Miller &; Scalo (1979) obtained much higher 
values of F = -0.4, -1.5, and -2.3 for, respectively, 0.1 < m < 1.0, 1.0 < m < 10, and m > 10. This 
suggests that R136's IMF is correlated with pi in a somewhat different way than the correlation 
that might exist in the Milky Way. 

In retrospect, this should not be surprising because R136 is much different than a spiral 
galaxy. The bright, early- type stars in spiral galaxies are generally confined to relatively narrow 
galactocentric radii near that of their initial birth sites. In contrast, stars in elliptical galaxies 
similar to R136 undergo substantial mixing due to their highly eccentric orbits. Therefore, the 
form of the IMF in spirals could be different than that in R136. 
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Fig. 1. — Top panel: the surface density of a spiral galaxy similar to the Milky Way but with the 
IMF of Model A. The dotted line is the surface density assuming that all stars lie on the main- 
sequence. The dashed line is the surface density if the V band mass to light ratio were constant 
at 7v = 2.0Mq/Lq. Lower three panels: circular velocities of Models A (upper-middle), B (lower- 
middle), and C (bottom). The circular velocities of Model A correspond to the surface density 
function shown in the top panel. For each model, the solid curve accounts for all components of 
mass, the dot-dash curve accounts for just the halo, dotted curve accounts for just the disk, the 
short-dashed curve accounts for just the bulge and spheroid stars, the dash-dotted curve accounts 
for everything except the halo, and the long-dashed curve assumes that the mass to light ratio is 
constant at tv = 2.0Mq/Lq. 
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For this reason, other models were also considered. Model B was constructed in order to help 

answer the question of just how necessary the dark halo is for circular velocity curves to be flat. The 
IMF gradient Fi was adjusted to minimize the curvature of the outer circular velocity curve, while 
Fq was adjusted such that the rotation velocity was ~ 220 km s~^. For simplicity, mi and ttIu were 
fixed. The middle plot of Figure 1 shows that the circular velocity curve of Model B is surprisingly 
fiat throughout most of the outer regions of the disk before the halo component is included. The 
total disk mass for Model B is 2.4x lO^^M©, with < 7v >disk =14:.5Mq/Lq, which is 7.3 times larger 
than the 7v = 2.0Mq/Lq model. The value of Fi for Model B is 0.42. This is 50% higher than the 
IMF gradient in R136. The change within the Milky Way of F measured by Garniany ct al. (1982) 
between the inner and outer semicircular regions of radius 2.5 kpc surrounding the Sun was -0.8, 
which for a disk scale length of i?o = 4.5 kpc corresponds to Fi = 0.8 x 37ri?o/(81ogex2.5 kpc)=3.9. 
This is much higher than the value in Model B. Thus the IMF gradient of Model B is well within 
empirical upper limits. 

However, there are at least five potential problems with the halo-less form of Model B: 1) In the 
adopted solar vicinity {Rq = 7.8 kpc), the surface density is 179M0pc~^. This is an unacceptable 
15 standard deviations higher than the local value of 46ib9 MqPC~^ measured by Kuijken Sz Gilmore 
(1989b). The corresponding mass to light ratio is 7.9Mq/Lq. This is 60% higher than the local 
value adopted in standard texts such as Binney &: Tremaine (1987). Similarly, the IMF slope at this 
radius is F = —1.7. Also, for the low mass (m ^ 0.5) stars, this value is incompatible with Miller & 
Scalo's (1979) result of F = —0.4. 2) Mestelian disks, which are similar to Model B, are commonly 
thought to be unstable to bar formation. The Toomre instability parameter Q is aiiK/{2.QGIline)j 
where k ~ 36 km s~^ kpc~^ is the epicycle frequency and is the mass-weighted stellar velocity 
dispersion (Toomre 1974) . Published estimates are Q ~ 1 — 3 in the solar vicinity. Because stellar 
velocity dispersions are empirically observed to decrease with mass even for stars with lifetimes 
greater than the age of the galaxy, estimates of aR are sensitive to mi. Wielen (1977) obtained 
aR = 62± 12 km s"^ for 0.1 ^ m ^ 0.8 K and M dwarfs, which implies Q ;^ 1.0 ± 0.2 for Model 
B. This lower limit is low enough to sustain spiral arm structure which numerical simulations show 
would rapidly dissipate otherwise. However, it is too near unity to prevent the growth of substantial 
arm/interarm stellar mass density contrasts. Though such mass contrasts are now known to exist 
in normal spirals (e.g., Rix Sz Zaritsky 1995, Gonzalez & Graham 1996), they are not accounted for 
in Model B. Incidentally, the halo component does not necessarily affect this instability (Sellwood 
1985). 3) The circular velocity curve at i? ;^ 35 kpc is not precisely flat, but actually rises before 
attaining a nearly Keplerian fall off. This is the result of the non-spherical potential. 4) The 
circular velocity drops below 200 km s~^ in the inner regions of the disk. This result is expected. 
For Model B, the mass to light ratio is ~ ml-y{mJmi)-'^-^{y + V){Mq/ Lq)/{-T - 1). If 

Fi = [ln(mu/mi) log^o e]~\ (1) 

which Model B obeys to within 15%, the mass to light ratio would scale as ~ This in turn 

would imply a disk surface density that is relatively constant. The circular velocities of such disks 
increase monotonically with R and are zero at R = 0. This problem with low inner disk velocities is 
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probably not serious because circular velocity curves are frequently compatible even with constant 
mass to light ratio, halo-less models throughout their entire optically-bright regions (e.g., Kent 
1986). Furthermore, flatter, halo- less velocity curves could probably be attained by including the 
following: galaxy parameters slightly different than those of Kuijken & Gilmore (1989a), a (more 
realistic) log-normal IMF (Miller & Scalo 1979), expected spatial dependence in remnant and gas 
mass fractions, and variations of mi or with pi. For instance, the velocity dip is better masked 
by the bulge if Bahcall & Soniera's (1984) smaller disk scale length of 3.5 kpc is assumed. 5) 
The IMF gradient of Model B appears to be too small to be compatible with the measurement 
of Garmany et al. (1982). Equation (|l|) suggests that this discrepancy would be less if a smaller 
m^/nii ratio had been employed. 

Of the above potential problems, only the first two appear to be significant at this time. Both 
can be overcome by taking into account the arm/interarm density contrasts observed in spiral 
galaxies; Model C was constructed to be similar to Model B, but has an azimuthally-averaged light 
and mass density that is 3.25 times greater than the interarm values in which the Sun presumably 
resides. The circular velocity curve of Model C, shown in the lower panel of Figure 1, is slightly 
higher but otherwise similar to that of Model B. However, the solar-vicinity disk surface density is 
only 6OM0PC-2. This is a much more reasonable 1.6 standard deviations above the value determined 
by Kuijken & Gilmore (1989b) and is actually lower than Bahcall & Soniera's (1984) value of 
~ 85M0PC-2. 

4. DISCUSSION 

A direct scaling of R136's IMF to the Galaxy does not dramatically alter the circular velocity 
curve. However, Models B and C, with their higher yet modest IMF gradients, have nearly flat 
"Wcirc ^ 220 km s~^ circular velocity curves only before the traditional dark halo component is 
included. It is interesting to note that if one assumes that these types of models and their ~ 10^- 
fold mass enhancements are representative of most galaxies, that the fiducial stellar contribution 
towards the closure density is ~ 0.004 (e.g., Peebles 1993) before accounting for IMF gradients, 
that the cosmological constant is zero, and that there is no hot dark matter, one would obtain 
Q ~ fibaryon — 0.04 -|- figas , where the closure fraction due to all gas including hot plasma in galactic 
clusters is 0.007 ^ ^gas^ 0.08 (Mulchaey et al. 1996). 

Current models of galactic evolution (e.g., Worthey 1994, de Jong 1996) do not account for 
IMFs that might vary with time and position via the temperature. This is despite prior warnings 
that the IMF probably has important dependences upon time and position (e.g., Mihalas &: Binney 
1978). In light of the above results, accounting for IMFs with such dependences may be necessary 
even to obtain results that are only accurate to first order. Accounting for these dependences 
may, for relatively obvious reasons, clarify our understanding of several astrophysical phenomena 
including the G-dwarf problem, intrinsic (as a function of radius) and extrinsic (as a function of 
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galactic morphology) metallicity and color gradients, and the Tully-Fisher relation. 

It is my pleasure to thank M. Harris, F. Bruhweiler, A. Whiting, A. Pridman, D. Kazanas, 
D. Audley, J. Knapen, A. Mignogna, E. Albert, and A. Stupp for useful comments. This work is 
part of a dissertation to be submitted to the Graduate School, University of Maryland in partial 
fulfillment of the requirements of the Ph.D. degree in Physics. 
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Table 1: IMFs in R136 



i?/pc 


r{R) 


m\ 


mu 


pi/(L0pc 3) 


0.20 


-1.29 ±0.20 


5.6 


120 


1.5 X 10^ 


0.60 


-1.46 ±0.23 


3.6 


76 


1.5 X 10^ 


2.0 


-2.12 ±0.09 


< 2.0 


48 


1.5 X 10^ 



Note. — Data adapted from Brandl et al. (1996) for (age-spread restricted) stars 2.5 — 3.5 Myr old. 



Table 2: Model Parameters 



Model 
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Ti 
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mil 
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Pme/ 




< Em© > / 


















M0PC-2 


Mopc-2 


A 


-3.03 


0.28 


-0.08 


0.13 


1.25 


0.12 


0.12 


67 


67 




(±0.26) 


(±0.06) 






(±0.23) 


(±0.04) 








B 


-1.11 


0.42 


-1.52 


0.00 


1.30 


0.00 


0.31 


179 


179 


C 


-0.55 


0.40 


-1.52 


0.00 


1.60 


0.00 


0.10 


60 


195 


Note. — 


This assumes f = fo 


+ /llogio[Pl/(-f'0PC" 


"')], for / 


= r, iogio('ni; 


), or logio(m 







